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ABSTRACT 


FURTHER  ONE-DIMENSIONAL  ANALYSIS  OF 
LONG-ROD  PENETRATION  OF 
SEMI-INFINITE  TARGETS 

by 

JOHN  DANIEL  CINNAMON,  B.S. 
SUPERVISING  PROFESSOR:  Wallace  T.  Fowler 


A  new  one-dimensional  analysis  of  long-rod 
penetration  of  semi-infinite  targets  is  presented.  Models 
in  this  field  attempt  to  accurately  describe  the  penetration 
process  of  a  long  rod  of  material  into  a  semi-infinite 
construct  of  target  material.  This  one-dimensional 
analysis  predicts  profile  hole  diameters  and  penetration 
depths  over  a  wide  range  of  material  combinations  and 
extends  the  previous  analysis  in  this  area  to  include 
hypervelocity  impacts.  This  approach  utilizes  values  of 
material  dynamic  yield  strengths,  known  impact 
conditions,  and  the  well  established  crater  volume/kinetic 
energy  relationship  to  predict  crater  hole  characteristics 
over  an  impact  velocity  range  of  1  to  6  km/s.  The  analysis 
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presented  here  includes  an  initial  transient  phase  and 
modifies  previous  estimates  for  pressure  on  the  penetrator 
tip  at  steady-state.  The  average  pressure  at  steady  state 
was  found  to  be  a  constant  value  over  the  range  of  impact 
velocities  for  a  particular  shot  combination.  The  specific 
value  for  the  average  pressure  was  found  to  be  a  direct 
function  of  target  strength.  The  resulting  equations  retain 
computational  simplicity  by  remaining  completely 
algebraic  in  nature.  Correlation  with  a  majority  of  the 
readily  available  experimental  data  is  given.  The  results 
are  quite  good  for  a  one-dimensional  model. 
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1.0  INTRODUCTION 


In  1967,  Tate  [1]  and  simultaneously  Alekseevskii  [2], 
published  a  one-dimensional  theory  for  the  penetration  of 
semi-infinite  targets  by  long  rods.  Tate  [3]  published  a 
second  paper  in  1969.  These  papers  form  the  basis  upon 
which  the  accepted  theory  of  one-dimensional  rod 
penetration  rests.  This  theory  attempts  to  accurately 
describe  the  penetration  process  of  a  long  rod  of  material 
into  a  semi-infinite  construct  of  target  material.  A 
prediction  of  the  resulting  penetration  hole  characteristics 
is  available  from  this  analysis.  By  comparing  the 
theoretical  predictions  to  experimental  data,  we  can 
evaluate  the  accuracy  of  these  one-dimensional  models. 

In  1987,  Jones,  et  al.  [4]  modified  the  equations  of 
motion  for  the  undeformed  rod  section  by  employing  a 
balance  of  linear  impulse  and  momentum.  In  subsequent 
papers  [5,6]  an  improvement  in  the  correlation  between 
theoretical  and  experimental  data  was  observed.  However, 
this  approach  relied  upon  post-test  measurements  for 
estimates  of  the  engineering  strain  in  the  mushroom  of  the 
penetrator  tip.  The  strain  was  assumed  to  be  constant 
throughout  the  process  and  equal  to  that  measured  from  the 
profile  hole  diameter.  By  using  the  modified  Bernoulli 
Equation  proposed  by  Tat6,  the  pressure,  penetration 
velocity,  and  undeformed  section  velocity  were  coupled 
with  the  modified  equation  of  motion  and  an  equation  for 
the  conservation  of  mass  passing  from  the  undeformed 
section  of  the  penetrator.  This  system  was  solved  by 
integration  and  the  resulting  prediction  of  penetration 
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depths  were  compared  to  experimentally  observed  values 
from  the  recovered  targets.  The  penetrator  and  target 
dynamic  yield  strengths  were  taken  to  be  constant  during 
the  event  and  estimated  from  laboratory  values  for  yield 
strengths  at  the  highest  available  strain-rates. 

These  results  were  satisfying  and  showed  promise, 
but  based  analysis  of  key  parameters  on  post-experiment 
measurement.  Kerber,  et  al.  [7],  in  1990,  utilized  a  well 
established  crater  volume  -  kinetic  energy  relationship, 
e.g.  [8],  to  remove  some  of  the  dependence  on  post-test 
measurements.  The  aim  was  to  produce  engineering  strain 
as  a  by-product  of  the  solution  of  the  model.  However, 
because  the  penetration  process  was  treated  as  being 
dominated  by  the  steady-state,  the  results  did  not 
correlate  well  over  a  large  range  of  impact  velocities  and 
rod  lengths.  The  correlations  with  experiment  were 
satisfactory  for  longer  rods  and  higher  impact  velocities 
where  steady  state  penetration  can  be  presumed  to 
dominate  the  event.  The  initial  and  terminal  transients 
were  neglected  in  this  and  previous  analysis.  The  initial 
transient  would  appear  to  become  more  important  in 
shorter  rods  and  lower  impact  velocities. 

In  1992,  Cinnamon,  et  al.  [9,10]  reported  a  significant 
increase  in  accuracy  by  incorporating  an  initial  transient 
phase  to  the  penetration  process.  This  analysis  was  based 
on  observations  by  Ravid,  et  al.  [11],  Gillis,  et  al.  [12],  and 
Jones,  et  al.  [13].  This  latest  approach  yielded  good 
accuracy  for  a  one-dimensional  model  in  the  1  to  3  km/s 
velocity  range.  The  results  depended  only  on  physical 
parameters  determined  before  testing  and  the  well 
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established  crater  volume/kinetic  energy  relationship.  The 
terminal  transient  was  neglected  in  this  model.  The 
resulting  equations  were  completely  algebraic  in  nature. 
Of  some  concern  in  this  approach  was  that  the  trends  in  the 
predicted  penetration  depth  curves  were  tending  toward 
significantly  high  values  for  impact  velocities  above 
3  km/s.  When  experimental  data  for  3  -  6  km/s  were 
evaluated,  the  model's  accuracy  indeed  deteriorated.  In 
examining  the  model,  it  became  obvious  that  the 

predictions  of  pressure  available  from  the  modified 

Bernoulli  Equation  were  simply  too  high.  In  fact,  the 

correlation  discovered  by  Cinnamon,  et  at.  [9,10]  indicated 
that  the  pressure  profile  could  be  successfully  modified  in 
such  a  way  that  the  average  pressure  at  steady-state 
predicted  by  the  Bernoulli  Equation  was  reduced  by  a 
specific  factor  (which  depended  exclusively  on  target 
strength)  to  match  experimental  penetration  depths.  The 
parabolic  nature  in  which  the  modified  Bernoulli  Equation 
predicts  the  interface  pressure  at  steady  state  as  impact 
velocity  increases  was  somewhat  minimized  by  this 

technique. 

In  1992,  Jones,  et  at.  [14],  the  pressure  distribution 
on  the  penetrator  tip  was  modified  to  attempt  to  correct 
for  the  deteriorating  results  above  2.5  km/s.  With  a 
judicious  choice  of  two  dimensionless  parameters,  slight 
improvement  was  observed.  However,  a  physical 
relationship  between  these  parameters  and  the  model  was 
not  discovered.  The  pressures  predicted  by  the  modified 
Bernoulli  Equation  were  once  again  reduced  to  achieve 
acceptable  results.  However,  this  new  pressure 


distribution  contained  the  same  parabolic  component  as  the 
previous  approach. 

In  this  paper,  several  approaches  to  the  choice  of 
pressure  at  steady  state  are  evaluated.  It  became  clear 
that  past  successful  solution  forms  for  the  pressure 
distribution  were  those  that  cancelled  out  the  parabolic 
nature  of  the  modified  Bernoulli  Equation.  As  a 
consequence,  when  the  average  pressure  at  steady  state  is 
taken  to  be  a  particular  value  for  a  specific  material 
combination  over  all  impact  velocities  and  directly  related 
to  target  strength,  a  significant  increase  in  model  accuracy 
resulted.  The  model  was  successfully  extended  up  to 
6  km/s.  This  paper  includes  a  large  body  of  results 
reflecting  analysis  of  all  readily  available  experimental 
data  in  this  field  of  research,  a  great  majority  of  which 
has  been  compiled  by  Anderson,  et  al.  [15]. 


2.0  THEORY 

The  general  concepts  of  the  rod  penetration  process 
are  detailed  in  Figure  A.  The  undeformed  penetrator  is  a 
cylinder  of  known  length  and  diameter  which  impacts  the 
target  at  a  nominal  normal  incidence  at  a  known  velocity. 
The  penetrator  enters  the  target  and  experiences 
mushrooming  in  the  tip.  When  the  event  has  concluded,  a 
crater  with  a  measurable  diameter  and  depth  remains. 
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Figure  A.  Schematic  of  Rod  and  Penetration  Process. 

(a)  shows  the  undeformed  rod  of  length  L  and  initial 
cross-sectional  area  Aj.  The  shaded  portion  will  be  lost 
to  erosion. 

(b)  shows  the  penetration  event.  I  is  the  undeformed 
section  length  and  z  is  the  penetration  depth. 

2.1  Primary  Governing  Equations 

In  [4],  Jones,  et  al.  proposed  a  modification  to  the 
equation  of  motion  of  the  undeformed  section  of  a  rod 
penetrator.  This  equation  is 


i;+i(v-u>  -  m 

where  l  is  the  undeformed  section  length,  v  is  the  current 

undeformed  section  velocity,  u  is  the  penetration  velocity, 
P  is  the  average  pressure  on  the  penetrator  tip,  p  is  the 
penetrator  density,  and  e  is  the  engineering  strain  in  the 
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penetrator  mushroom.  The  details  behind  the  development 
of  this  equation  are  contained  in  [4]. 

Wilson,  et  al.  in  [5]  added  another  key  equation  to  the 
analysis.  The  conservation  of  mass  across  the  plastic 
interface  between  the  mushroom  and  the  undeformed 
section  of  the  penetrator  was  given  in  the  form 


et  =  v-u  (2) 

The  penetrator  is  assumed  to  be  rigid-plastic  during  the 
event.  In  both  (1)  and  (2),  dots  over  the  symbols  represent 
differentiation  with  respect  to  time.  The  engineering 
strain  in  the  mushroom  is  compressive  and  therefore 
negative.  The  current  value  of  the  engineering  strain  in  the 
mushroom  is  defined  to  be 

As 

e  =  ^  -  1  (3) 

where  Aj  is  the  initial  cross-sectional  area  of  the 
undeformed  penetrator  and  A  is  the  current  cross-sectional 
area  of  the  mushroom. 

2.2  Previous  Pressure  Analysis 

The  average  pressure  on  the  penetrator  tip,  P,  can  be 
varied  by  considering  various  pressure  profiles,  in  [4]  and 
[5],  the  pressure  was  assumed  to  be  uniform  across  the 
mushroom  face.  The  intensity  of  the  pressure  was  assumed 
to  be  the  solution  of  the  modified  Bernoulli  equation.  This 
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equation,  from  [1  -  3],  is  applied  at  steady  state  and 
relates  pressure  on  the  axis  of  the  specimen  at  the 
penetrator  tip,  pa,  to  the  undeformed  section  speed  v,  the 

penetration  velocity  u,  and  material  properties  of  the 
target  and  penetrator.  The  modified  Bernoulli  equation  is 

Paa52M2pu2  +  Rt=  2  p  (V'u>2  +  Yp  (4) 

where  Rt  and  Yp  are  dynamic  yield  strengths  of  the  target 
and  penetrator  at  suitably  high  strain  rates  respectively, 
and  m  2  is  the  ratio  of  the  penetrator  to  target  density.  The 
pressure  was  taken  in  [4]  and  [5]  not  to  vary  across  the 
mushroom  face. 

Gillis,  et  al.  ,  in  [12]  suggested  that  the  uniform 
pressure  profile  was  not  realistic  and  proposed  a  parabolic 
form  which  was  symmetric  about  the  axis  of  the  specimen 
and  zero  on  the  edge  of  the  mushroom.  In  this  case,  the 
pressure  p  had  the  form 


P  -  Pa  (1  -  )  (5) 

where  pa  is  the  current  pressure  on  the  rod  axis,  R  is  the 
original  rod  radius,  and  r  is  the  radial  distance  from  the 
axis.  The  factor  (1+e)  in  the  denominator  of  the  pressure 
term  in  (1)  forces  the  pressure  to  act  over  the  deformed 
mushroom  face  with  area  A,  even  though  (5)  refers  to  the 
original  rod  configuration. 

The  average  pressure  P  can  be  computed,  in  general, 
by 


8 


P 


J  P  dAj 

As 


(6) 


Pa 

When  P  was  calculated  for  (5),  the  result  was  ^  ■  The 

predictions  for  penetration  depths  were  somewhat 
improved  over  the  previously  assumed  uniform  pressure 
*  distribution.  The  results  still  suffered  from  ignoring  the 
initial  and  final  transients  in  the  penetration  event.  This 
parabolic  pressure  distribution  effectively  reduced  the 
parabolic  nature  of  (4)  by  a  factor  of  two.  Although  this 
improved  the  model  by  reducing  the  effect  of  (4),  the 
addition  of  at  least  an  initial  transient  seemed  warranted. 

In  Cinnamon,  et  al.  [9,10],  the  pressure  distribution 
was  generalized  to  the  form 


P 


Pa  (1 


kn 


R2 


(7) 


which  made  the  average  pressure  term  become 


P  = 


_Pa_ 

(n+1) 


(8) 


This  new  pressure  distribution  was  successfully 
employed  with  an  initial  transient  phase  to  correlate  to  a 
large  number  of  experimental  cases.  The  pressure  exponent 
n  was  found  to  be  a  direct  function  of  target  strength. 
Although  this  approach  provided^  improved  results  over 
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previous  approaches,  the  analysis  could  noi  be  successfully 
extended  beyond  3  km/s. 

Jones,  et  al.  [14]  attempted  to  correct  the  problem  by 
adding  a  uniform  component  to  the  pressure  profile.  The 
distribution  had  the  form 

p  =  q  +  (pa  -  q)  (1  -  )n  (9) 


and  the  average  pressure  became 


(10) 


This  new  pressure  profile  improved  results  slightly  and 
extended  them  into  the  hypervelocity  range  (3-6  km/s). 
However,  q  and  n  were  not  successfully  correlated  to  any 
physical  parameters. 

In  this  paper,  another  approach  to  the  determination 
of  the  average  pressure  will  be  explored  and  a  successful 
correlation  to  target  strength  will  be  presented.  In 
addition,  the  widest  possible  body  of  available  data  will  be 
used  in  the  correlation. 


2.3  Transient  Penetration  Analysis 

In  Cinnamon,  et  al.  [9,10],  an  initial  transient  phase 
of  penetration  was  added  to  the  model.  The  transient  phase 
is  characterized  by  impact  shock  effects  and  complete 
mushroom  growth  which  occurs  between  impact  and  the 
beginning  of  steady  state  penetration.  We  assume  that  the 
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penetrator  impacts  the  target  at  a  known  velocity,  v0,  of 
sufficient  magnitude  (i.e.  v0  >  1  km/s)  such  that  the 
undeformed  section  cannot  sustain  any  appreciable 

deceleration  (i.e.  v  s  0)  during  the  initial  transient.  This 
means  that  v  *  v0  throughout  the  transient  or  mushrooming 

phase  of  penetration.  During  the  initial  transient,  the 
mushroom  develops  from  a  cross-sectional  area  A0  at 
impact  to  A1t  when  steady  state  begins.  The  mushroom 
retains  an  area  of  A1  throughout  the  steady  state  portion  of 
the  penetration  process  until  the  end  of  the  event. 

Ravid,  et  al.  [11]  reported  that  there  was  little 
change  in  penetration  velocity  u  during  the  shock/impact 
stage  of  the  initial  transient  phase.  Motivated  by  this 
observation,  we  assumed  that  the  penetration  velocity  was 
approximately  constant  (i.e.  u  =  u0)  throughout  the  initial 

transient.  Hence,  (1)  becomes 


(ID 

and  (2)  is  modified  to 

m 

el  =  v0  -  u0  (12) 

These  two  equations,  (11)  and  (12),  govern  the 
mushrooming  of  the  rod  during  the  initial  transient  phase 
of  the  penetration  event  which  precedes  the  steady  state. 
At  impact,  the  engineering  strain  in  the  mushroom  is  e0. 
When  steady  state  is  reached,  the  strain  becomes  e^ 
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By  eliminating  l  between  (11)  and  (12)  and  solving 
for  e,  we  arrive  at  an  expression  for  the  engineering  strain. 


-(Vq-Uq)2 
(V0“U0)2  +  “ 


(13) 


This  relationship  provides  us  with  an  explicit  formula  for 
the  strain  in  the  mushroom  as  it  develops  during  the 
transient  phase.  This  equation  governs  the  behavior  until 
the  beginning  of  the  steady  state  penetration  phase.  The 
pressure  on  the  axis,  pa,  is  changing  rapidly  during 
mushroom  formation.  It  has  a  large  value,  p0,  at  impact 
and  a  reduced  value,  pt ,  at  steady  state. 

When  steady  state  is  reached,  we  assume  that  the 
modified  Bernoulli  Equation,  (4),  is  valid.  At  the  transition 
point  between  the  transient  and  steady  state  portions  of 
the  event,  (4)  can  be  expressed  as 

Pi  =  £M2P  u02  +  R,=  l  p  (v0  -  u0)2  +  Yp  (14) 


and  (13)  can  be  written  as 


-(Vq-Uq)2 

(v0-u0)2  +  ^ 


(15) 


where  P1  is  the  average  pressure  on  the  penetrator  tip  at 
the  beginning  of  steady  state. 
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Equation  (14)  can  then  be  used  to  solve  for  u0  in 
terms  of  the  known  quantities  v0,  p,  jm2,  Rt,  and  Yp.  The 
remaining  variable  in  the  system  of  equations  is  Pj.  The 
determination  of  P-j  is  the  primary  focus  of  this  paper. 

The  penetration  velocity  u0  can  be  found  algebraically 
from  (14).  The  three  primary  cases  are  outlined  below. 

For  equal  penetrator  and  target  dynamic  yield 
strength  (i.e.  Rt=Yp)  and  equal  densities  (i.e.  m2=1),  u0 

reduces  to 


u0»^v0  (16) 

For  unequal  penetrator  and  target  dynamic  yield 
strength  (i.e.  Yp)  and  equal  densities  (i.e.  m2*1),  u0 

becomes 


uo 


pvn2+2(Yp-R,) 

2pv0 


In  the  general  case,  u0  is  given  by 


(17) 


-V, 


1/2 


(18) 


In  each  of  these  cases,  then,  u0  can  be  determined 

algebraically  from  known  material  properties  and  impact 
conditions.  This  leaves  e-|  in  (15)  as  a  function  of  known 

parameters  and  the  average  pressure  on  the  penetrator  tip 
at  steady  state. 
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2.4  Impact  Conditions 

The  model  outlined  above  also  provides  some 
information  about  conditions  at  impact.  The  strain  on 
impact  e0  can  be  calculated  from  (13)  if  we  know  the 
average  pressure  on  the  penetrator  tip  at  impact,  P0. 


-(Vq-Uq)2 

(v0-u0)2  +  ^ 


(19) 


The  impact  pressure  can  be  estimated  from  elementary 
shock  physics,  using 


Po  -  P  us  u0  (20) 

where  us  is  the  shock  speed  in  the  target.  Values  for  us  as 
a  function  of  u0  can  be  found  in  shock  Hugoniot  tables,  e.g. 
[16].  Calculation  of  P0  from  a  known  value  of  p0  can 

typically  be  accomplished  using  the  same  approach  as  the 
calculation  of  P j  from  pa. 


2.5  Cratering  Analysis 

The  mathematical  model  for  the  behavior  of  the 
penetrator  is  a  rigid-plastic,  instantaneously  eroding  rod 
model.  As  a  result,  the  penetrator  enters  the  target  with 
some  impact  engineering  strain  e0  which  expands  to  e1 
during  the  transient.  The  impact  pressure  p0  is  usually 
very  high  relative  to  the  steady  state  pressure  p1 . 
Although  this  pressure  decreases  rapidly  during  mushroom 
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formation  in  the  transient  phase,  the  values  for  p0  can  be 
significant.  The  mushroom  diameter  grows  from  the  time 
of  impact  through  the  transient  phase,  and  ceases  at  the 
beginning  of  the  steady  state  portion  of  the  event.  The 
shock/impact  stage  takes  place  in  a  period  of  a  few 
microseconds,  see  [11]. 

The  instantaneous  erosion  assumption  prevents  the 
model  from  accounting  for  any  additional  erosion  of  the 
target  -  which  occurs  in  actual  practice.  There  is  typically 
appreciable  change  in  target  geometry  due  to  penetrator 
and  target  material  ejection  from  the  crater.  As  a 
consequence,  the  recovered  targets  will  appear  to  have 
more  cylindrical-type  craters  than  the  model  would 
predict.  Figure  B  illustrates  the  crater  predicted  by  the 
mathematical  model,  and  Figure  C  indicates  how  the  actual 
geometry  frequently  appears. 


Figure  B.  Idealized  Crater  Geometry 

A0  and  A1  are  the  cross-sectional  areas  at  impact  and  at 
steady  state  respectively,  z  is  the  penetration  depth. 
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Figure  C.  Actual  Crater  Geometry 

A1  is  the  observed  crater  cross-sectional  area.  The 
crater  is  assumed  cylindrical  with  altitude  z. 


2.6  Penetration  Analysis 

Predicting  penetration  depths  from  the  above  model 
is  made  possible  through  the  use  of  a  somewhat  empirical 
approach.  For  a  number  of  years,  researchers  have  observed 
a  significant  correlation  between  crater  volume  in  the 
recovered  targets  and  impact  kinetic  energy,  e.g.  [8].  This 
relationship  appears  to  be  linear  for  impact  cases  of 
sufficiently,  but  not  excessively,  high  energy  and  can  be 
expressed  in  the  form 


Vc  =  a  E0  +  b  (21) 

where  Vc  is  the  crater  volume,  E0  is  the  impact  kinetic 
energy  and  is  equal  to  ^pAjLv02,  and  the  variables  a  and  b 
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are  regression  constants  determined  from  the  available 
experimental  data.  The  linear  fit  is  performed  for  each 
shot  combination. 

These  crater  volurne/kinetic  energy  relationships  are 
computed  from  data  points  for  a  particular  shot 
combination.  The  reliability  of  the  linear  fit  is,  of  course, 
a  function  of  the  number  of  data  points  available.  Since 
most  experimental  tests  are  quite  expensive,  frequently 
the  data  is  sparse  and/or  somewhat  scattered.  When  the 
linear  fit  predicts  cratering  at  zero  impact  kinetic  energy, 
or  in  some  other  way  reflects  erroneous  trends,  the 
usefulness  of  the  particular  case  is  significantly  reduced. 
This  phenomenon  is  typically  avoided  by  a  sufficient 
number  of  experimental  points.  The  accuracy  of  the 
penetration  prediction  is  extremely  dependent  on  the  crater 
volume/kinetic  energy  relationship  arrived  at  using  this 
technique. 

By  adopting  the  cylindrical  approximation  for  the 
crater  geometry  discussed  above,  we  can  generate 
predictions  for  the  penetration  depths.  Because  of  the 
ejection  of  material  from  the  crater,  the  cross-sectional 
area  of  the  recovered  target  hole,  AlP  will  be 


A-, 

Ai  ■  (i  +  ®i > 
The  crater  volume  can  be  expressed  as 


(22) 


Vc  -  A,  z 


(23) 


where  z  is  the  penetration  depth.  The  penetration  depth 
can  be  predicted  by  applying  the  crater  volume/kinetic 
energy  relationship  to  Vc.  From  (21)  and  (23),  z  is  given  by 


z  =  f  (1+e,)  Vc  =  ^-(1  +ei)(aE0  +  b)  (24) 

Thus,  the  penetration  depth  can  be  expressed  as  an 
algebraic  function  of  known  material  properties  and  impact 
conditions,  and  as  a  function  of  the  average  pressure  on  the 
penetrator  tip  at  steady  state. 


3.0  CURRENT  PRESSURE  PROFILE  ANALYSIS 

In  the  previous  work  outlined  above,  it  was  noted  that 
one  of  the  primary  difficulties  in  achieving  good 
predictions  for  the  crater  characteristics  was  the  manner 
in  which  the  modified  Bernoulli  Equation  (4)  predicted 
pressure  as  a  function  of  u  and  v.  The  parabolic  nature  of 
(4)  tended  to  over-predict  penetration  depths  for  the  higher 
velocity  cases  (3-6  km/s).  Earlier  successes, 
particularly  in  the  intermediate  velocity  range  (1  - 
3  km/s),  were  the  result  of  pressure  distributions  that 
tended  to  reduce  the  effect  of  (4).  In  this  paper,  additional 
pressure  distribution  analysis  is  performed  to  address  this 
problem. 
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3.1  Previous  Pressure  Distributions 

An  obvious  departure  point  for  the  attempt  to 
improve  the  one-dimensional  analysis  and  extend  it  into 
the  hypervelocity  range  was  to  begin  with  the  previous 
pressure  distributions.  A  great  body  of  additional 
experimental  data  became  available  in  [15]  that  expanded 
the  range  of  materials,  greatly  increased  the  number  of 
shot  combinations  accessible  for  analysis,  and  provided 
data  in  the  hypervelocity  range.  The  pressure  profiles 
detailed  in  (7)  and  (9)  were  examined  for  possible 
application  in  these  new  cases. 

The  distribution  in  (7)  was  unable  to  compensate  for 
the  parabolic  nature  of  (4)  while  maintaining  accuracy  in 
the  intermediate  velocity  ranges  if  n  was  considered  a 
constant  as  it  was  in  [9,10].  Attempts  to  model  n  as  a 
function  of  material  properties  or  impact  conditions  were 
not  successful.  In  general,  it  was  observed  that  n  needed 
to  increase  with  impact  velocity  to  essentially  cancel  the 
effects  of  the  modified  Bernoulli  Equation.  The 
distribution  in  (9)  followed  the  same  trend.  Although  the 
results  improved,  the  net  effect  was  to  choose  n  and  q  to 
counter  the  dramatic  pressure  increase  dictated  by 
equation  (4). 

3.2  New  Pressure  Distribution 

To  attempt  to  find  a  solution  to  this  dilemna,  a  new, 
more  general,  pressure  distribution  was  proposed  by  Jones 
[17].  This  pressure  profile  takes  the  form 


This  profile  is  a  more  complex  and  versatile  one.  A  great 
deal  of  additional  control  over  the  shape  of  the  pressure 
distribution  was  provided  by  (25).  The  average  pressure 
then  is  given  by 

r<m)  r(1+n) 

P  =  q  +  2<pa-q) - j -  (26) 

m  r(1+im+n) 

where  r  is  the  well  known  mathematical  gamma  function. 

As  a  significant  number  of  cases  were  examined,  it 
became  increasingly  clear  that  in  order  to  achieve  the 
desired  trends  in  the  theoretical  penetration  curves,  n  and 
m  were  chosen  in  such  a  way  as  to  essentially  eliminate 
the  effect  of  the  second  term  in  (26).  That  is,  an  average 
pressure  comprised  of  a  single  value,  q,  which  was 
unvarying  over  the  range  of  impact  velocities,  achieved  the 
best  results.  This  discovery  matched  our  previous 
experience  with  (7)  and  (9).  Apparently,  the  magnitude  and 
trends  in  the  pressures  predicted  by  the  modified  Bernoulli 
Equation  were  not  leading  to  acceptable  results.  All 
previous  successes  were  based  on  choices  that  reduced  or 
eliminated  the  contribution  of  the  velocity  dependent  axial 
pressure  in  (4)  to  the  value  of 


3.3  Revised  Average  Pressure  Approach 

With  the  results  described  above,  another  approach 
was  required.  It  became  clear  that  the  previous  calculation 
of  the  average  pressure  at  steady  state,  P-j,  was  not 
acceptable.  To  simplify  the  analysis,  the  connection  of  P-, 
to  a  particular  pressure  distribution  is  ignored. 

The  problem  is  further  simplified  by  the  fact  that  the 
desired  trends  in  the  penetration  depths  result  from  a 
constant  value  for  P1  over  the  entire  velocity  range. 

This  approach  does  not  imply  that  the  pressure 
distributions  have  the  same  shape  for  differing  velocities, 
or  that  pa  is  equivalent  for  all  velocities,  simply  that  the 
value  of  P1  is  constant  for  a  particular  shot  combination 
over  all  impact  velocities. 


4.0  RESULTS 

When  P1  was  assumed  to  be  some  constant  average 
pressure  on  the  penetrator  tip  independent  of  v0,  defined 
here  to  be  Q,  for  all  impact  velocities  in  a  particular  shot 
combination,  the  results  of  this  model  improved 
tremendously.  The  penetration  depth  theoretical  curves 
adopted  the  trends  present  in  the  experimental  data  (i.e. 
penetration  depths  leveling  off  as  velocity  increases 
toward  6  km/s).  In  addition,  the  values  for  Q  which  yielded 
the  best  results  correlated  strongly  to  target  strength. 

In  order  to  establish  the  most  credible  and  most 
complete  correlation  possible,  all  available  data  were 
employed.  As  a  consequence,  this  report  includes  a  large 
number  of  figures.  In  order  for  this  theory  to  be  applied 
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effectively,  experimental  data  sets  must  have  included 
crater  diameters.  A  number  of  the  cases  in  [15]  did  not 
provide  this  information.  In  addition,  a  minimum  of  two 
data  points  were  required  to  construct  the  crater- 
volume/kinetic  energy  relationship.  Hence,  some  other 
cases  could  not  be  evaluated.  With  those  limitations  in 
mind,  the  author  applied  all  readily  available  cases  to  this 
model  and  reports  the  results. 

Table  1  summarizes  all  the  cases  and  provides 
essential  data  about  each  shot  combination.  The  figure 
numbers  referred  to  can  be  found  in  Appendix  A.  One  figure 
number  is  assigned  to  each  shot  combination.  For  each 
case,  several  figures  denoted  by  the  figure  number  and 
lower  case  letters  are  to  be  found.  The  first  three  figures 
of  each  case  (corresponding  to  a,  b,  and  c)  are  the  strain  vs. 
impact  velocity,  crater  volume  vs.  impact  kinetic  energy, 
and  penetration  depth  vs.  impact  velocity  graphs.  Any 
additional  graphs  will  occur  in  pairs.  These  additional 
figures  (i.e.  d  and  e,  etc.)  refer  to  cases  in  which  the  crater 
volume/kinetic  energy  relationship  was  modified  by 
removing  certain  data  points  that  appear  erroneous  or  that 
may  skew  the  general  trend.  When  these  points  are 
removed,  the  relationship  changes  -  which  modifies  the 
resulting  penetration  curve.  Of  course  this  does  not  alter 
the  strain  curve  in  any  way  -  avoiding  the  necessity  of 
including  an  additional  one.  The  first  three  figures  for 
each  case  include  all  reported  data. 

The  figures  report  theoretical  curves  superimposed 
on  discrete  experimental  points.  The  upper  curve  in  the 
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Table  1.  Experimental  Data  Summary 


Fig. 

# 

Pene- 

trator 

YP 

(MPa) 

P 

1%) 

L 

(mm) 

L 

D 

Target 

Rt 

(MPa) 

«2e 

m° 

Ref. 

1 

Alum 

Alloy 

200 

2700 

63.5 

10 

Lead 

200 

11200 

[18] 

2 

1100- 
0  Al 

250 

2720 

9.525 

3 

1100- 
0  Al 

250 

2720 

[19] 

3 

2024- 
T3  Al 

675 

2770 

9.525 

3 

1100- 
O  Al 

250 

2720 

[19] 

4 

7075- 
T6  Al 

600 

2804 

9.525 

3 

1100- 
O  Al 

250 

2720 

[19] 

5 

C1015 

St. 

600 

7600 

9.525 

3 

1100- 
O  Al 

250 

2720 

[19] 

6 

Soft 

4340 

1263 

7850 

31.75 

5 

2024- 
T4  Al 

400 

2770 

[20] 

7 

7075- 
T6  Al 

600 

2804 

31.75 

5 

2024- 
T4  Al 

400 

2770 

[20] 

8 

Hard 

4340 

1826 

7850 

31.75 

5 

7075- 
T6  Al 

600 

2804 

[20] 

9 

Soft 

4340 

1263 

7850 

31.75 

5 

7075- 
T6  Al 

600 

2804 

[20] 

10 

7075- 
T6  Al 

600 

2804 

31.75 

5 

7075- 
T6  Al 

600 

2804 

[20] 

11 

1100- 
0  Al 

250 

2720 

9.525 

3 

Cl  01 5 
St. 

600 

7600 

[19] 

12 

2024- 
T3  Al 

675 

2770 

9.525 

3 

C1015 

St. 

600 

7600 

[19] 

13 

7075- 
T6  Al 

600 

2804 

9.525 

3 

C1015 

St. 

600 

7600 

[19] 

14 

C1015 

St. 

600 

7600 

9.525 

3 

C1015 

St. 

600 

7600 

[19] 

15 

4340 

Steel 

1600 

7810 

68.58 

18 

6061- 

T651 

Al 

600 

2710 

[21] 

18 

1  100- 
0  Al 

250 

2720 

9.525 

3 

2024- 
T3  Al 

675 

2770 

[19] 

Pena-  Yp 
trator 


(MPa)  kg. .  (mm) 

l  <i ) 


(MPa) 


17 

1100- 
O  Al 

250 

2720 

9.525 

3 

304 

St.St. 

675 

18 

2024- 
T3  Al 

675 

2770 

9.525 

3 

2024- 
T3  Al 

675 

19 

2024- 
T3  Al 

675 

2770 

i 

9.525 

3 

304 

St.St. 

675 

20 

m 

600 

7600 

9.525 

3 

2024- 
T3  Al 

675 

21 

C1015 

St. 

600 

7600 

9.525 

3 

KfjQ| 

Tma 

675 

22 

304 

St.St. 

675 

7900 

9.525 

3 

2024- 
T3  Al 

675 

23 

304 

St.St. 

675 

7900 

9.525 

3 

304 
St.  St. 

675 

24 

C110 

W1 

1200 

7850 

25 

10 

St37 

750 

25 

C110 

W1 

1200 

7850 

43 

10 

St37 

750 

26 

C110 

W1 

1200 

7850 

54 

1 

1 

St37 

750 

27 

C110 

W1 

1200 

7850 

25 

10 

St52 

850 

28 

C110 

W1 

1200 

7850 

43 

10 

St52 

850 

29 

Clio 

W1 

1200 

7850 

54 

10 

St52 

850 

30 

D17 

2500 

in 

28 

10 

St52 

850 

31 

D17 

2500 

IKZUlQI 

60 

10 

St52 

850 

32 

C110 

W2 

m 

7850 

58 

10 

HzB.A 

850 

33 

Marag 

St 

1000 

7850 

58 

1 

| 

HzB.A 

850 

34 

Marag 

St 

1000 

7850 

116 

20 

HzB.A 

850 

35 

35CrN 

iMo 

2200 

7850 

54 

10 

HzB.A 

850 

36 

elmet 

MWMMU 

15500 

58 

10 

HzB.A 

850 

7900 


2770 


7900 


2770 


7900 


2770 


7900 


785C 


7850 


7850 


7850 


7850 


7850 


1 1531 

\vm 


ii^i 
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Fig. 

# 

Pene- 

trator 

YP 

(MPa) 

M 

m 

L 

(mm) 

| 

Target 

Rt 

(MPa) 

m2p 

m*3 

Ref. 

SH 

fPffl 

37 

D17K 

2500 

17300 

58 

wm 

HzB.A 

850 

7850 

vm 

38 

l>IM 

2500 

17000 

58 

10 

HzB.A 

850 

7850 

fESl 

39 

W 

2500 

19300 

60 

10 

HzB.A 

850 

7850 

Wm 

40 

W75 

2500 

15500 

58 

10 

HzB.A 

850 

7850 

vm 

41 

W90 

2500 

17000 

58 

10 

HzB.A 

850 

7850 

pm 

42 

D18 

2500 

18000 

58 

10 

HzB.A 

850 

7850 

1031 

sxa 

H01T 

2000 

14500 

58 

10 

HzB.A 

850 

7850 

wm 

i  44 

D17.6 

C&Ittl 

1 7600 

41.7 

10 

HzB.A 

850 

7850 

wm 

warn 

D17.6 

EM 

17600 

58 

10 

HzB.A 

850 

7850 

IES1 

KOI 

Steel 

1600 

7850 

58 

10 

HzBA 

850 

7850 

wm 

wm 

D17 

2500 

17000 

116 

W4*  ■ 

HzB.A 

850 
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strain  vs.  impact  velocity  figures  represents  the  estimate 
for  e0. 

When  these  cases  were  evaluated,  a  certain  value  for 
Q  could  be  chosen  to  match  the  experimental  data.  This  Q 
was  found  to  strongly  correlate  to  target  strength.  Table  2 
reports  each  of  the  different  targets  present  in  the  data 
and  their  corresponding  Q  value. 
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Table  2.  Correlation  of  Q  to  Target  Strength 


Target 

Rt 

Q 

Material 

(MPa) 

(GPa) 

Lead 

200 

0.128 

1100-0  Al 

250 

0.28 

2024-T4  Al 

400 

0.6 

7075-T6  Al 

600 

1.05 

Cl  01 5  St. 

600 

1.05 

6061-T651  Al 

600 

1.05 

2024-T3  Al 

675 

1.3 

304  St.  St. 

675 

1.3 

St37 

750 

1.4 

St52 

850 

1.7 

HzB.A 

850 

1.7 

RHA 

1000 

2.2 

W8 

1000 

2.2 

Ger  RHA 

1000 

2.2 

Ger  Arm  St 

1100 

2.275 

Ger  St 

1200 

2.35 

Soft  4340 

1263 

2.4 

4340  Steel 

1600 

2.5 

Hard  4340 

1826 

2.6 

D17 

2500 

2.9 

With  this  information,  it  was  immediately  evident 
that  some  strong  relationship  existed  between  Q  and  target 
strength.  Figure  D  depicts  the  values  for  Q  chosen  to  allow 
the  model  to  predict  penetration  depths  and  crater 
diameters  accurately  against  Rt.  A  curve  is  fit  through  the 

data  to  both  illustrate  the  correlation  and  provide  a 
functional  relationship  between  Q  and  Rt.  The  best  fit  is 


Q  «  3.8  (1  -  e  (*0  00135  Rt)  )  -  0.8 


(27) 
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0-3  00135*Rt))-  8 


Target  Strength  (GPo) 


figure  D.  Average  Pressure  vs.  Target  Strength 
Q  is  the  average  pressure. 


With  Q  as  a  direct  function  of  target  strength,  this 
one-dimensional  penetration  model  can  be  expressed  in 
terms  of  known  material  properties  and  impact  conditions. 
The  crater  volume/kinetic  energy  curve  is  still  needed  to 
allow  for  the  calculation  of  penetration  depths,  however. 

To  place  this  revised  model  into  context  with  regard 
to  those  used  in  [9,10]  and  [14],  Figures  E  and  F  are 
presented.  These  figures  illustrate  the  comparative 
results  between  these  three  different  approaches  to  the 
calculation  of  the  average  pressure  at  steady  state. 


Impact  Velocity  (km/*) 

Figure  E.  Comparative  Case:  Strain  vs.  Impact  Velocity 

(a)  indicates  theory  from  [9,10] 

(b)  indicates  theory  from  [14] 

(c)  indicates  current  model 


Kennemetol  WIO  on  RHA 


hnpoct  velocity  (km/e) 

Figure  F.  Comparative  Case:  Penetration  Depth  vs.  impact  Velocity 

(a)  indicates  theory  from  [9,10] 

(b)  inefcates  theory  from  [14] 

(c)  indicates  current  model 


It  is  clear  from  the  above  comparative  figures  (i.e.  E 
&  F)  that  each  subsequent  theory  lowered  the  average 
pressure  in  such  a  way  as  to  lower  the  penetration  and 
strain  curves.  This  particular  case  that  appears  in  Figures 
E  and  F  was  chosen  for  its  higher  velocity  data  and  its 
presentation  in  [14]. 


5.0  CONCLUSION 

In  this  paper,  a  new  one-dimensional  model  for  the 
penetration  of  semi-infinite  targets  by  long  rods  was 
developed.  Its  basis  was  the  revision  of  the  previous 
techniques  employed  to  calculate  the  average  pressure  at 
steady  state.  The  modified  Bernoulli  Equation  appears  to 
result  in  values  for  that  are  simply  too  high  for  the 

impact  cases  in  the  hypervelocity  range. 

By  correlating  a  new  approach  to  a  great  body  of  data, 
it  was  discovered  that  a  single  value  for  the  average 
pressure  on  the  penetrator  tip  at  steady  state,  Q,  could 
successfully  represent  the  pressure  at  steady  state  over 
the  impact  velocity  range  of  1  to  6  km/s.  This  formulation 
for  Q  allowed  the  model  to  improve  its  accuracy  and  its 
trends  at  higher  velocities.  In  addition,  this  value  Q  was 
shown  to  have  a  strong  correlation  to  target  strength.  This 
approach  has  resulted  in  a  completely  algebraic  solution 
that  relies  only  on  known  test  parameters  and  the  well 
established  crater  volume/kinetic  energy  relationship. 
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Future  work  will  involve  an  effort  to  revise  or 
replace  the  modified  Bernoulli  Equation's  estimate  for 
pressure  at  steady  state.  Additional  analysis  also  needs  to 
be  conducted  to  ascertain  the  form  of  the  pressure 
distribution  that  leads  to  a  constant  Q  over  all  impact 
velocities. 

The  aim  of  this  paper  was  to  extend  the  one¬ 
dimensional  penetration  analysis  into  the  hypervelocity 
range.  The  resulting  model  offers  good  accuracy  for  a  one¬ 
dimensional  description  of  the  penetration  event. 


Appendix  A: 
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Figure  la.  Strain  vs.  Impact  Velocity 
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Figure  1b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  1c.  Penetration  Depth  vs.  Impact  Velocity 


AJ  urn  in  urn  AJtoy  on  L*od 


Vc  »1. 282 33*Eo— 1405 


•mpact  Kinetic  Energy  (IcJ) 

Figure  Id.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  2b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  2c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  5a.  Strain  vs.  Impact  Velocity 
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Figure  5b.  Crater  Volume  vs.  impact  Kinetic  Energy 
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Figure  6b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  6c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  7a.  Strain  vs.  Impact  Velocity 
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Figure  7b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  7c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  7d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  8b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  8c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  8d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  8e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  9e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  10a.  Strain  vs.  Impact  Velocity 
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Figure  11a.  Strain  vs.  Impact  Velocity 
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Figure  11b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  11c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  12a.  Strain  vs.  Impact  Velocity 
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Figure  15a.  Strain  vs.  Impact  Velocity 
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Figure  15b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  15c.  Penetration  Depth  vs.  Impact 
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Figure  16a.  Strain  vs.  Impact  Velocity 
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Figure  16b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  16c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  17a.  Strain  vs.  Impact  Velocity 
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Figure  17b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  18b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  18c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  19a.  Strain  vs.  Impact  Velocity 
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Figure  19b.  Crater  Volume  vs.  Impact  Kinetic  Energy 


impact  Velocity  (km/«) 


Impact  Kinetic  Energy  (kj) 


Impact  Velocity  (km/*) 


Figure  21a.  Strain  vs.  Impact  Velocity 


C101S  Steel  on  JO*  Stoimei*  Steel 


Figure  21b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  22b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  22c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  25a.  Strain  vs.  Impact  Velocity 
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Figure  25b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  25c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  25d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  26b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  2$c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  27a.  Strain  vs.  Impact  Velocity 
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Figure  27b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  27e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  28a.  Strain  vs.  Impact  Velocity 
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Figure  28b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  2|c.  Penetration  Depth  vs.  Impact  Velocity 


•mpoct  Velocity  (km/ s) 

Figure  29a.  Strain  vs.  Impact  Velocity 
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Figure  29b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  30b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  30c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  31c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  32b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  32c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  32d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  32e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  34b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  34c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  35b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  35c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  35d.  Crater  Volume  vs.  impact  Kinetic  Energy 
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Figure  35e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  36a.  Strain  vs.  Impact  Velocity 
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Figure  36b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  36c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  37c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  37e.  Penetration  Depth  vs.  Impact 
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Figure  38a.  Strain  vs.  Impact  Velocity 
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Figure  38b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  38c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  38d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  38e.  Penetration  Depth  vs.  Impact  Velocity 


(cpuosrtoijx) 

(£%.ujuj)  •ujnjoA 


104 


0 

-0.1 

-0.2 
-0.3 
-0.* 

-0  5 
-0  6 
-0.7 
-0.8 
-0.9 

-1 

0  2  4  6 

impact  Velocity  (km/i) 

Figure  39a.  Strain  vs.  Impact  Velocity 
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Figure  39b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  39c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  39d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  41a.  Strain  vs.  Impact  Velocity 
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Figure  42b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  42d.  Crater  Volume  vs.  Impact  Kinetic  Energy 


D18  on  HtBA 


import  Velocity  (km/*) 


Figure  42e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  44b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  44d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  45a.  Strain  vs.  Impact  Velocity 
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Figure  45b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  45c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  45d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  45e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  45f.  Crater  Volume  vs.  Impact  Kinetic  Energy 


Impact  Klnatlc  Enargy  (kJ) 


Figure  46b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  47a.  Strain  vs.  Impact  Velocity 
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Figure  47b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  48b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  48c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  49a.  Strain  vs.  Impact  Velocity 
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Figure  49b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  49c.  Penetration  Depth  vs.  Impact  Velocity 


Figure  49d.  Crater  Volume  vs.  Impact  Kinetic  Energy 


Impact  Velocity  (km/») 


Figure  49e.  Penetration  Depth  vs.  Impact  Velocity 


H60T  on  H2B.A 


impoct  Velocity  (km/a) 

Figure  50a.  Strain  vs.  Impact  Velocity 


Impact  VaiocKy  (km/*) 


Figure  51a.  Strain  vs.  Impact  Velocity 
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Figure  52b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  53b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  55a.  Strain  vs.  Impact  Velocity 
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Figure  55b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  57a.  Strain  vs.  Impact  Velocity 
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Figure  59a.  Strain  vs.  Impact  Velocity 
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Figure  59b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  61b.  Crater  Volume  vs.  Impact  Kinetic  Energy 


Strain  P#n*t ration  Oapth  (mm) 


110 
100 
90 
60 
70 
60 
50 
40 
30 
20 
10 
0 

0  2  4  6 

Impact  Velocity  (km/#) 

Figure  61c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  62a.  Strain  vs.  Impact  Velocity 
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Figure  62b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  62d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  62e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  63c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  64c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  65a.  Strain  vs.  Impact  Velocity 
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Figure  65c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  66d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  67a.  Strain  vs.  Impact  Velocity 
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Figure  67b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  68b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  68c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  68d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  68e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  69b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  70b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  70d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  71c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  73b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  73c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  73d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  75b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  75e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  76a.  Strain  vs.  Impact  Velocity 
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Figure  76b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  76c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  77a.  Strain  vs.  Impact  Velocity 
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Figure  77b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  77c.  Penetration  Depth  vs.  impact  Velocity 
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Figure  77e.  Penetration  Depth  vs.  Impact  Velocity 


0 

-0.1 

-0.2 
-0.3 
-0.4 
-0.5 
-0  6 
-0.7 
-0.8 

-0.9 

-1 

0  2  4  6 

Impact  Velocity  (km/e) 

Figure  78a.  Strain  vs.  Impact  Velocity 
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Figure  79a.  Strain  vs.  Impact  Velocity 
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Figure  79b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  80d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  80e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  81a.  Strain  vs.  Impact  Velocity 
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Figure  81b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  82b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  82d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  82e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  83c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  83d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  84a.  Strain  vs.  Impact  Velocity 
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Figure  84b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  84d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  85a.  Strain  vs.  Impact  Velocity 
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Figure  85b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  85e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  86a.  Strain  vs.  Impact  Velocity 
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Figure  86b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  88b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  88c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  88d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  88e.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  89a.  Strain  vs.  Impact  Velocity 
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Figure  89b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  90b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  90c.  Penetration  Depth  vs.  Impact  Velocity 
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Figure  91a.  Strain  vs.  impact  Velocity 
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Figure  91b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  91  d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  93a.  Strain  vs.  Impact  Velocity 
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Figure  93b.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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Figure  94d.  Crater  Volume  vs.  Impact  Kinetic  Energy 
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